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ABSTRACT 

Tij" \ We report on the analysis of the first complete far- infrared spectrum (43-197/^m) of the Seyfert 2 

O^l ■ galaxy NGC 1068 as observed with the Long Wavelength Spectrometer (LWS) onboard the Infrared 

' Space Observatory (ISO). In addition to the 7 expected ionic fine structure emission lines, the OH 

\ rotational lines at 79, 119 and 163/xm were all detected in emission, which is unique among galaxies 

, with full LWS spectra, where the 119/xm line, when detected, is always in absorption. The observed 

line intensities were modelled together with ISO Short Wavelength Spectrometer (SWS) and optical 
and ultraviolet line intensities from the literature, considering two independent emission components: 
the AGN component and the starburst component in the circumnuclear ring of 3 kpc in size. 
C^h' Using the UV to mid-IR emission line spectrum to constrain the nuclear ionizing continuum, we 

have confirmed previous results: a canonical power-law ionizing spectrum is a poorer fit than one 
with a deep absorption trough, while the presence of a big blue bump is ruled out. Based on the 
^ instantaneous starburst age of 5 Myr constrained by the Br 7 equivalent width in the starburst 

\ ring, and starburst synthesis models of the mid- and far-infrared fine-structure line emission, a low 

ionization parameter (U=10~^'^) and low densities (n=100 cm~'^) are derived. Combining the AGN 
and starburst components, we succeeded in modeling the overall UV to far-IR atomic spectrum of 
NGC 1068, reproducing the line fluxes to within a factor 2.0 on average with a standard deviation 
of 1.3, and the overall continuum as the sum of the contribution of the thermal dust emission in 
' the ionized and neutral components. The OH 119 /im emission indicates that the line is coUisionally 

excited, and arises in a warm and dense region. The OH emission has been modeled using spherically 
symmetric, non-local, non-LTE radiative transfer models. The models indicate that the bulk of the 
emission arises from the nuclear region, although some extended contribution from the starburst is 
not ruled out. The OH abundance in the nuclear region is expected to be ^ 10~^, characteristic of 
X-ray dominated regions. 

Subject headings: galaxies: individual (NGC 1068) - galaxies: active - galaxies: nuclei - galaxies: 
Seyfert - galaxies: emission lines - galaxies: starburst - infrared: galaxies. 

1. INTRODUCTION 

NGC 1068 is known as the archetypical Seyfert type 2 galaxy. It is nearby, luminous (Lir ~ 2 x IO^^Lq 
iBland- Hawthorn et al.lir997D . and it has been e xtensively o bserved and studied in detail from X-rays to radio wave- 
lengths. With a measured redshift of z=0.0038 l)Huchra et al. 1999.) (corresponding to a distance of D=15.2 Mpc for 
Hn=75 km Mp c~^), it provides a scale of only ~ 74 pc/". A cei itral nuclear star cluster has an extent of 0.6" 
ijThatte et al.ll997j) and a 2.3 kpc stellar b ar observed in the near- IR (IScoville et aUllQS^lThronson et al.llT989t) IS sur- 
rounded by a circumnuclear starburst ring. .Telesco et alJ l|T984) found that the infrared emission in NGC 1068 was due 
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to both the Seyfert nucleus (which dominates the 10/xm emission) and to the star forming regi ons in the bright ~ 3kpc 
circumnuclear ring (which emits most of the luminosity at A > 30^m). A Br^ imaging study (jPavies. Sugai fc Ward I 
Il998j) showed a similar morphology and indicated that a short burst of star formation occurred throughout the cir- 
cumnuclear ring of 15-16" in radius within the last 4-40 Myr. CO interferometer observation s revealed molecular ga s 
very close to the nucleus (^0.2") suggesting the presence of ~ 10® within the central 25pc ('Schinn erer et al.l l2000'). 
Recent high resolution H2 line emission mapping indicates the presence of two main nuclear emission knots with a 
velocity difference of 140 km/s, which, if interpreted as quasi-keplerian, would imply a central enclosed mass of IQ^Mq 
l|Alloin et al.ll200l . 

In this article, we present the first complete far-infrared spectrum from 43 to 197/im showing both atomic and 
molecular emission lines (§2). We model the composite UV- to far-IR atomic emission line and continuum spectrum, 
from our data and the literature, using photoionization models of both the active nucleus and the starburst component 
(§3). We also model the mid- to far-IR continuum emission using a radiative transfer code and gray body functions 
for the neutral molecular components (§4). Moreover, two different non-local, non-LTE radiative transfer codes have 
been used to model the OH lines (§5). Our conclusions are then given in §6. 

2. OBSERVATIONS 

NGC 1068 was observed with the Long Wav elength Spectrometer (LWS) l|Clegg et alJ ll996"l on board the Infrared 
Space Observatory (ISO) l)Kessler et aljTl99^ . as part of the Guaranteed Time Programme of the LWS instrument 
team. The full low resolution spectrum (43-197/j,m) of NGC 1068 was collected during orbit 605 (July 13, 1997). Two 
on-source full scans (15,730 seconds of total integration time) and two off-source (6' N) scans of the [CII]158/im line 
(3,390 seconds of total integration time) were obtained. On- and off-source scans had the same integration time per 
spectral step. Because of the design of the LWS spectrometer, simultaneously with the 158^m data, a short spectral 
scan of equal sensitivity to the on-source spectrum was obtained at sparsely spaced wavelengths across the LWS range. 

The LWS beam is roughly independent of wavelength and equ a.1 to about 80 arcsec. The spectra were calibrated 
using Uranus, resulting in an absolute accuracy better than 30% (jSwinvard et al. 1996'). The data analysis has been 
done with ISAP^, starting from the auto-analysis results processed through the LWS Version 7-8 pipeline (July 1998). 
To be confident that newer versions of the pipeline and calibration files did not yield different results, we have compared 
our data with the results obtained using pipeline 10.1 (November 2001) and we did not find significant differences in 
the line fluxes or the continuum. 

All the full grating scans taken on the on-source position and the two sets of data on the off-source position were 
separately co-added. No signal was detected in the off-source coadds. The emission line fluxes were measured with 
ISAP, which fits polynomials to the local continuum and Gaussian profiles to the lines. In all cases the observed line 
widths were consistent with the instrumental resolution of the grating, which was typically 1500 km/sec. The integrated 
line fluxes measured independently from data taken in the two scan directions agreed very well, to within 10%. The on 
source LWS spectrum that resulted from stitching the ten LWS channels together using small multiplicative corrections 
in order to match the overlapping regions of each channel with its neighbors is shown in Fig. 1. LWS spectra of sources 
that are very extended wit hin the instrument bea m or that peak off center are typically affected by channel fringing 
in the continuum baseline ijSwinvard et al.| 11998(1 . Fortunately, these spurious ripples are hardly noticeable in our 
LWS spectrum, presumably because the far-IR continuum is centrally concentrated towards the center of the LWS 80" 
beam. 

■ ! . ' . :• the LWS observations, we also use the SWS observations presented bv iLutz et al.l i200r)D . to extend the 
wavelength and ionization-level coverage. Table ^ presents all the ISO line flux measurements including those from 
the SWS with their respective aperture sizes. 

3. THE FINE STRUCTURE LINES 

To be able to better constrain the modeling of the line emission of NGC 1068, we have combined our far-infrared 
fine structure line measurements (Tabled) with ultraviolet, optical and infr ared spectroscopic data from the literature 
l|Kriss et al.lll992HMarconi et al.ltT996t iThomnson II1996I: iLutz et al.ll200fD . The complete emission line spectrum of 
NGC 1068 from the ultraviolet to the far-IR includes several low-ionization lines that are primarily produced outside 
the narrow line region (NLR) of the active nucleus, as well as intermediate ionization lines that originate from both 
starburst and AGN emission. For this reason, we find that no single model satisfactorily explains all the observed 
emission lines. We identify two main components: 

- an AGN component (the NLR), exciting the high ionization lines and contributing little to the low-to-intermediate 
ionization lines; 

- a starburst component in the circumnuclear ring of the galaxy fe.g. lDavies. Sugai fc Ward 1 11998') that produces 
the low ionization and neutral forbidden lines and some of the emission in the intermediate ioniz ation lines. This 
comp onent should also produce emission associated with photo-dissociation regions (PDRs) fe.g. lKaufman et alJ 
II999D . at the interface with the interstellar medium of the galaxy. 

^ The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument Teams and Data Centers. 
Contributing institutes are Centre d'Etudc Spatiale des Rayonnements (France), Institute d'Astrophysique Spatiale (France), Infrared 
Processing and Analysis Center (United States), Max-Planck- Insitut fiir Extraterrestrische Physisk (Germany), Rutherford Apploton 
Laboratories United Kingdom) and the Space Research Organization, Netherlands. 
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In this section, we will examine separately the two components that produce the total fine structure emission line 
spectrum of NGC 1068, namely the AGN and the starburst, for which we propose two different computations, and we 
add together these components to reproduce the overall observed spectrum from the UV to the far-lR in §3.3. 

3.1. Modeling the AGN 

The first photoionzation model predictions of the mid to far-infrared emission line spectra of the Narrow Line Regions 
(NLR) of active galaxies w e re pre sented by[Spino2lio & Malkan (1992), well before the ISO observations could be 
collected. [Alexander et alJ (f20W) used the observed high ionization emission lines to model the obscured ionizing 
AGN continuum of NGC 1068 and found that the best-fit spectral energy distribution (SED) has a deep trough at 
4 Rydbergs, which is consistent with an intrinsic "big blue bump" that is partially obscured by ^ 6 x 10^^ cm^^ 
of neutral hydrogen interior to the NLR. Following their results, we have simulated th eir models, although using a 
different photoionization code, CLOUDY {Yersion 96. lFerland et alJll998HFerland II2000D and t hen we have varied the 
shape of the ionizing continuum to include the ionizing continuum derived in iPier et all |1994). Our goal was to test 
if the Alexander et al,: (2000) results were unique and to fit the remaining emission by a starburst component, and 
thereby to derive a composite model of the complete emission line spectrum of NGC 1068. 

Specifically, w e explo re three plausible AGN SEDs. Model A assumes the best fit ionizing spectrum derived by 
[Alexander et all (|2000t). i.e. with a deep trough at 4 Rydberg (log/ = -27.4, -29.0, -27. 4, -28.2 at 2, 4, 8 and 16 
Ryd, respectively). An intrinsic nuclear spectrum of NGC 1068 has also been inferred bv iPier et al.l fl994'). Model 
B assumes the original io nizing spectrum d erived from Pier et al. ( 1994,) . Model C assumes an SED with a Big Blue 
Bump superposed on the IPier et all ({19941 ^ ionizing continuum (log/ = -25.8, -25.8, -25.8, -27.4 at 2, 4, 8 and 16 
Ryd, respectively) as expected for the thermal emission of an accretion disk around a central black hole. These three 
AGN ionizing continua are plot ted in Fig. U For each of models A, B, and C, we have used two component models 
with the same parameters as in [Alexander et akl ()2000() : component 1 has a constant hydrogen density of lO^'cm^'^, 
an ionization parameter U=0.1, a covering factor c=0.45, a filling factor of 6.5 x 10^^ with a radial dependence of 
the form r~^, and extends from ^ 21 to 109 pc from the center; component 2 has a density of 2 x 10'^ cm~^, an 
ionization parameter U=0.01, a covering factor of c=0.29, a filling factor of 6.5 x 10"'* without any radial dependence, 
and extends from ^ 15 3 to ~ 362 pc from the center. We have also assumed the "low oxygen" abundances adopted by 
[Alexander et al.l 1)20001) . in order to be able to compare our results with theirs^. Because the grain physics has been 
updated in the most recent version of CLOUDY (Version 96), we have included the presence of grains in the models, 
using "Orion- type" grains'^. The inclusion of grains also allows us to compute the thermal dust continuum emission 
from the ionized components (see §4). 

The inner and outer radii of the emission regions of the two components, 21, 109, 153 and 362 pc, correspond 
to angular distances of about 0.26, 1.4, 1.9 and 4.5 ", respectively. Table 12 reports the predicted line fluxes of the 
three AGN models. A, B, and C, together with the observed line fluxes: the line fluxes are given for each of the two 
components 1 and 2, which are treated as independent, and the total flux for each model is simply the sum of the 
fluxes of the two components. 

We can see from Table |21 that only the AGN A and B models, and not the AGN C model, reproduce most of the 
observed high ionization line fluxes. The low and intermediate ionization lines, are expected to have partial or full 
contributions from starburst and PDR components (sec §3.2). This first result rules out the presence of a "big blue 
bump" in the ionizing continuum of NGC 1068. To be able to compare the modeled ultraviolet and optical lines with 
the observations, we al so listed in Tab le El and 01 their dereddened fl uxes, assumi ng two values for the extinction: 
Es-y = 0.4 mag (Mal kan fc Oke II1983D and ^b-v = 0.2 mag (Marcon i et al. Ill 9961) . We find that the AGN B model 
overpredicts several of the intermediate ionization lines, such as [SIV]10.5/im, [NeIII]15.6/im and [SIII]18.7/im, and 
this discrepancy increases when adding the starburst component because these lines are also copiously produced by 
that component (see next section). On the other hand, the [NeII]12.8/im emission is underprcdicted so much so that 
even with the inclusion of the starburst component it cannot be reproduced with this model. As we discuss further in 
§3.3, a composite AGN/starburst model using AGN model A reproduces the [NeII]12.8/im emission better than the 
other composite models. 

3.2. Modeling the starburst ring 

NGC 1068 is known to emit strong starburst emission from the ring-like structure a t a radial distance of 15 — 16" 
from the nucleus (total size of 3 kpc), traced for example by the Br 7 emission (iDavies. Sueai fc Ward 1 11998). 
Mid-IR line imaging observations of NGC 1068 have been published by iLe Floc'h et all l|200lD based on~OCAM 
CVF observations. They presented an image of the 7.7 /im PAH feature that shows constant surface brightness above 
the 4th coutour near the nucleus. This suggests that star formation is occuring in the direction of the nucleus so 
that nuclear spectra will include some emis sion from star formation. In the case of the SWS observations that we 
are modeling (reported bv lLutz et alJ l|2000D 'l. three apertures were used at different wavelengths with the two largest 
also including portions of the brighter starburst ring (see Table To estimate how much of the starburst emission 

^ The adopted gas phase chemical abundances in logarithmic form are: H: 0.00, He: -1.00, Li: -8.69, Be: -10.58, B: -9.21, C:-3.43, N:-3.96, 
O : -3.57, F : -7.52, Ne: -3.96, Na: -5.67, Mg: -4.43, Al: -5.53, Si: -4.46, P : -6.49, S : -4.79, CI: -6.72, Ar: -5.60, K : -6.88, Ca: -5.64, Sc: 
-8.83, Ti: -6.98, V : -8.00, Cr: -6.33, Mn: -6.54, Fe: -4.40, Co: -7.08, Ni: -5.75, Cu: -7.79, Zn: -7.40 

^ The abundances of the grain chemical composition, in logarithmic form, are: C: -3.6259, O: -3.9526, Mg: -4.5547, Si: -4.5547, Fe: 
-4.5547 
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is contained in the different apertures used in the observations, we have used a continuum subtracted image in the 
6.2 Mm feature p roduced by C. Dudley (private communication) using the same ISOCAM CVF data set examined by 
ILe Floc'h et"alfll200lh . The 6.2 fj.-m feature is more isolated than the 7.7 ^m feature, which is blended with the 8.7 /xm 
feature and the silicate absorption feature, but the image compares well with the published 7.7 fim image though we 
have zeroed out the residuals in a 3x12 arcsec^ region centered on the nucleus. Based on this image, the SWS 14x20, 
14x27 and 20x33 arcsec^ slits contain 13, 23 and 46% of the 6.2 fim flux contained in the LWS beam respectively, 
without correction for the neglected region of poor residuals (oriented at 45° to our synthetic SWS slits). Since PAH 
features are thought to be a good tracer of PDRs and their associated startbursts, we adopt these percentages in our 
model predictions of SWS line strengths in the starburst models presented in this section. In fitting our starburst 
models to the observations, we have computed the line fluxes at earth of each centrally illuminated emitting cloud and 
then determined the number of clouds needed to best fit the observed li ne fluxes. 

We have chosen the starburst synthesis modeling program Starburst99 ijLeitherer et al.ll99flt) . to produce input ioniz- 
ing spectral energy distributions (SEDs) for the CLOC/£>Fphotoionization code. Models were followed to temperatures 
down to 50 K to include the PDR components. We compared the predictions of an instantaneous star formation law 
with those of a continuous star formation law. For both types of models we adopted an age of 5 Myr, a Salpeter IMF 
(a=2.35), a lower cut-off mass of 1 Mq, an upper cut-off mass of 100 Mq, solar abundances (Z=0.020) and nebular 
emission included. These ionizing SEDs are shown in Fig. |21 with a total mass of M = 10^ Mq for the instantaneous 
model and a star formation rate of 1 Mq yr~^ for the continuous model. These partic ular ioni zing continuum shapes 
were selected because they are consistent with the Br 7 equivalent width observed bv lDavie~ Sugai & Ward (199^ 
in the starburst ring. We have estimated that the Br 7 equivalent width in each of the individual regions of the map 
of Da vies. Sugai fc Ward (1998 ) is in the range 110-180 A. According to the Lcithcrcr et al. (1999) models (see their 
figures 89 and 90), for a value of log(W(Br 7, A)) > 2 only instantaneous models with ages less than ~ 6x 10^ yrs 
are allowed. 

We report in Table 13 the line fluxes predicted for six different centrally illuminated starburst models, chosing the 
above instantaneous star formation SED as the input ionizing continuum and using CLOUDY with densities of iih = 
10, 100, 1000 cm~^, ionization parameters of Log U = -2.5, -3.5 and an inner cloud radius of 50 pc. As a function of the 
adopted density, we then determined the following numbers of emitting clouds needed to fit the observations: 33000, 
3300, and 330 clouds for the three values of the density, respectively. The fluxes reported in Table 3 are therefore 
the total starburst line fluxes at earth and together with the nuclear line fluxes of Table 2 can be compared with the 
observations. 

We have also run models with the continuous star formation law presented above, but we do not list their results in 
Tabled because the differences in the line flux predictions compared with the instantaneous models are insignificant, 
compared with the effects of density and ionization parameter, as can be seen from TableOl This result is not surprising 
because the ionizing continua of the two starburst models are quite similar in shape and the total number of clouds 
is a free parameter. We have also tried continuous starburst models with much longer ages (10, 20 and 100 x 10^ 
years) but, because the shape of the ionizing continuum again does not change significantly, the resulting emission line 
spectrum was indistinguishable from that one derived from the models with an age of 5 x 10^ years. 

In all models the abundances were those typical of HII regions^ and grains of "Orion-type"^ are included. The 
integration was allowed to run until the temperature of the gas in the cloud cooled to T=50 K in order to include the 
photodissociation regions present at the interfaces of HII regions and molecular clouds. 

It is clear from Table 01 that the models with the higher ionization parameter (log U — -2.5) can easily be ruled 
out, because their emission in many intermediate ionization lines is far too high (see e.g. [0IV]26/J,m, [OIII]51,88/im, 
[NIII]57/im). Among the models with the lower ionization parameter (log U = -3.5), we can exclude model SBR F, 
with density — 1000 cm"'^, because it underestimates many far-IR lines which are not strongly emitted by the 
active nucleus (namely: [SiII]35/im, [NIII]57/im, [OIII]88/im, [NII]122/im, [CII]158/im) while the low density model 
(SBR B, with uh — 10 cm^'^) overpredicts the [CII]158/im line by a factor of 2 relative to the other far-IR lines and 
does not reproduce the [OIII] doublet ratio. Finally, the intermediate density model (SBR D, with uh = 100 cm~^) 
gives the best fit to the observed lines, taking into account that the AGN component must be added to reproduce the 
total flux as is shown in §3.3. 

Wc estimate the average PDR parameters using the models of iKaufman et all ()1999D and the contour plots in 
iLuhman et al. (200^, the measured [C II]158 and [O I]145 ^m fine fluxes (but not the [O I]63 fim line flux which 
may be affected by absorption and/or shocks), and the FIR flux integrated over the LWS spectrum, which we find to 
be 1.3 X 10~^ ergs cm~^ sec~^. Here we assume that the [C II] line emerges predominantly from PDRs due to the 
strong starburst, rather than the diffuse ionized medium. With this assumption, the average PDR gas density and UV 
radiation field are njj^ ~ 1000 and Go ~ 300 respectively. We note that if instead we assume that the [C II] line flux is 
dominated by the diffuse ionized medium, using the correction factor estimated by,Malhotra ct al. (2001J), we obtain 
a similar gas density njjj ^ 1500 but a significantly higher interstellar radiation field Gn ~ 1500. For both cases, the 
parameters derived are in the range of those of the normal galaxies in the .Malhotra et al.l l|200lD sample, consistent 

The adopted gas phase chemical abundances in logarithmic form are: H: 0.00, He: -1.02, Li: -10.27, Be: -20.00, B: -10.05, C: -3.52, N: 
-4.15, O: -3.40, F: -20.00, Ne: -4.22, Na: -6.52, Mg: -5.52, Al: -6.70, Si: -5.40, P: -6.80, S: -5.00, CI: -7.00, Ar: -5.52, K: -7.96, Ca: -7.70, 
Sc:-20.00, Ti: -9.24, V:-10.00, Cr: -8.00, Mn: -7.64, Fe: -5.52, Co:-20.00, Ni: -7.00, Cu: -8.82, Zn: -7.70 

^ The abundances of the grain chemical composition in logarithmic form are: C : -3.3249, O : -3.6516, Mg: -4.2537, Si: -4.2537, Fe: 
-4.2537. 
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with the assumption that most of the FIR flux originates in the starburst ring. 

3.3. Adding the two components 

Summing the hue intensities of each one of the two components, the composite spectrum of NGC 1068 can be 
derived and compared with the observed one. We have chosen three combinations to compute the composite models, 
each one with a different AGN model, while we adopted the starbust model with nn = 100 cm~'^ and Log U = -3.5: 
1) the first one (that we name CMl, for Composite Model 1) with the AG N ionizing c ontinuum as suggested by 
lAlexander et all ^QOd) (model AGN A); 2) the second (CM2) with the original ' Pie'r et all l|1994) (model AGN B); 3) 
the third (CM3) with the hypothetical bump (model AGN C). The results of these three composite models are given in 
Table0I compared to the observed and dereddened values, assuming the two choices for the extinction (see §3.1). We 
also show the results of the three composite models in a graphical way in Fig. |31 where the modeled to the observed 
flux ratio is given for each line for the case of an extinction of Eb-v=0.2 mag. 

A simple test of the three models resulted in a reduced of 11.6, 17.1 and 177 for the three models CMl, CM2 
and CM3, respectively, for an extinction of Fib-v—OA mag, while these values become 23, 46 and 325 for EB_y=0.2 
mag. Thus, of the models explored, CMl with Eb_\/=0.4 mag provides our best flt to the observations. We note that 
model CMl reproduces the line fluxes to within a factor of 2.0 on average, with a standard deviation of 1.3. 

4. THERMAL CONTINUUM SPECTRUM 

To model the total mid- and far-infrared thermal dust continuum of NGC 1068 we have used different computations 
for each component. The modelled emission is shown for the individual components and combined in Fig. 01 The 
thermal dust emission from the AGN narrow line regions and the starburst regions in the ring have been computed using 
the CLOUDY photoionization models described in the previous sections. Specifically, the UV continuum reprocessed 
by the dust present in both NLR components 1 and 2 of model AGN A has been diluted by the same covering factors 
that affect the line emission (c=0.45 and 0.29 for components 1 and 2, respectively). Similarly we computed the 
continuum emission from our best starburst model (SBR D). While wc find that the continuum produced in this way 
for the AGN is consistent with the observed mid-infrared energy distribution (Fig. ^J, accounting for about half of 
the observed emission at mid-infrared wavelengths, the emission from dust associated with the starburst ionized and 
photodissociated regions, although similar in shape to the observed continuum, produces only ~ 20% of the far-IR 
continuum. We have indeed performed a search in parameter space by varying the age of the starburst (from 4 to 6 
X 10^ yrs), the gas density (from 10 to 1000 cm~^) and the radius of the emitting clouds (from 25 to 100 pc), but 
no starburst model that could reproduce simultaneously the observed line and far-infrared continuum emission was 
found. 

If the CLOUDY models correctly refiect conditions in both the ionized and photodissociated gas, these results 
may imply that the bulk of the starburst far-infrared continuum arises from dust that is mixed with neutral gas 
not directly associated with the ionized or photodissociated gas. However, because photoionization codes such as 
CLOUDY have not been used in the past to model the dust continuum, and because we may not have fully searched 
parameter space, we are hesitant to over-interpret these results until more detailed comparison with galactic ionized 
and photoionized regions are car ried out. We have therefore described this starburst thermal dust component, following 
ISpinoglio. Andreani fc Malkan 1 12002), in terms of a gray body function with a temperature of T = 34 K, and a colder 
gray body component at T = 20 K (see Fig. 0J. These are gray body functions with a steep (/3 = 2) dust emissivity 
law. Assuming a spherical shell with radius of 1.5 kpc and thickness of 0.3 kpc, the inferred average H2 molecular 
density associated with the 34 K component is 3.7 cm~^. The total mass is 1.2 x 10^ Mq and 2 x 10^ Mq for the 34 
K and 20 K components, respectively. T hese estimates are in reasonable agreement with the 4 x 10^ Mq derived 
for the molecular ring from CO emission l|Planesas et alJll991|) . 

As pointed out above, the dust associated with the ionized NLR components 1 & 2 of model AGN A does not quite 
account for the total mid-infrared continuum (Fig. 0J. Therefore, we have assumed that the missin g mid-infrared 
arises from the ne utral-nuclear component, which has been observed in a variety of molecular lines (e.g iTacconi et all 
Il994t l y sero et aTl 12004). For this neutral component we have assumed a total gas+dust mass of M = 2 x 10^ 
M^ l|Helfer fc Bl itz ||199 5[l. and modelled the expected con tinuum using a non-local, spherically symmetric, radiative 
transfer code jConzdlez^lfonso fc Cernicharo ||19 97lll999f) . The molecular nuclear emission has been resolved into a 
circumnuclear disk or ring l)Schimere^^^l2000() . which here is roughly modelled as a dusty spherical envelope with 
inner and outer radii of 3 and 200 pc, respectively. We assume an AGN luminosity of 3.7 x 10^" Lq, which is a factor 
of ~ 3 lower than the total AGN luminosity, simulating that most of the AGN radiation escapes thro ugh the poles of 
the m olecular disk, and/or is absorbed in the NLR, and so is not able to heat the molecular gas (e.g. iCameron et all 
Il993j) . The dust envelope is divided into a set of spherical shells where the dust temperature is computed assuming 
that heating and cooling are equal. We used a standard silicate/amorphous carbon mixture with optical constants 
given bv lDrain'el l)1985|) and lPreibisch et all l)1993|) . The density profile was assumed to be oc r~^ , with /? regarded as 
a free parameter. 

The resulting mid-infrared emission, obtained with /3 = 1, is shown in Fig. 01 The nuclear molecular component has 
an averaged H2 density at the inner radius of < 717(112) >= 500 cm~^, and a column density of A^(H2) = 2 x 10^^ cm~^. 
Once the emission from this neutral component is summed up with the emission predicted for the NLR components 
1 & 2, a good fit to the observed emission for A > 9 fim is obtained. For A < 9 /im, the mid-infrared emission is 
underestimated, suggesting the presence of a hot component, probably very close to the central AGN, which is not 
included in our models. 
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5. THE OH LINES 
5.1. General remarks 

In NGC 1068, we detect three of the OH rotational lines, all in emission. As shown in the energy level diagram of 
Fig-El two of them are fundamental lines, connecting the ground state ^113/23/2 level with the ^113/25/2 (the in-ladder 
119 fj,m line) and with the ^ni/2l/2 level (the cross ladder 79 fim line). The third line is the lowest transition of the 
ladder: the 163/im line between the J=3/2 and J=l/2 levels. The detected line fluxes are given in Tabled 
The fact that these three lines are all in emission is in striking contra st with t he OH lines observed i n other bright 
infrared g alaxies, such as Arp 2 20 (jPischer et al. 1999; Gonzalez- Alfons o et 311 12004) . Mrk 231 (Harv ev et alJ ll9991. 
NGC 253 l|Bradford et al.lll999D . and M 82 itColbert et alJll999 1. in which the 119 ^ni fundamental is in absorption. 
The 79/im line is sometimes seen in emission and sometimes in absorption; the 163/.tm line is always seen in emission. 
In addition to the detections, the ISO-LWS and SWS observations provide upper limits on fluxes of the other four 
lines that arise between the lowest six rotational levels. The LWS spectra in the vicinity of the detected lines (and of 
one of the upper limits), are shown in detail in Fig. |^ (histograms). In this section we discuss the physical conditions 
necessary to excite these lines, their probable location within NGC 1068, and detailed model fits to their fiuxes. A 
comparison between the observed and modeled line fluxes is given in Table |5l and shown in Fig.|H| 

5.2. The excitation mechanism of the OH lines 

The unique OH emission line spectrum of NGC 1068 can provide a powerful way to help discriminate between the 
properties of the molecular clouds in NGC 1068 and the clouds in other galaxies in which OH has been observed. 
Before describing our detailed radiative transfer calculations, it is instructive to discuss some conclusions that are 
model- independent. The emission in the OH n3/2 5/2 — 3/2 line at 119 /im cannot be explained by absorption of 
far- infrared photons followed by cascade down to the upper H3/25/2 level of the transition. Rather, we argue that 
collisional excitation dominates. Figure [S] shows the energy level diagram of OH. There are only two possible paths 
to excite the 119 /im line via absorption of far-infrared photons: via the 35 and/or the 53 ground-state lines. 
Excitation by either of these routes has other observable consequences. In the case of simple radiative cascading, the 
Einstein-A coefficients of the lines involved in the corresponding cascades are such that if the 35 /im absorption path 
were responsible for the observed 119/im line flux, then the OH H1/2 5/2 — 3/2 line at 98.7 /im would be approximately 
5 times stronger than the 119 /im line, while the 98.7 /im line is not detected. Hence this possibility is ruled out. 
Similarly, if absorption in the 53/im line were responsible for the observed 119/im line flux, then the 163/im line would 
be about 5 times stronger than the 119 /im line, which it is not. We can therefore conclude from the constraints 
provided by the other far-infrared OH lines that the 119/im emission line is not the result of radiative absorption 
and cascading. The implication is that OH excitation through collisions is more important in NGC 1068 than in the 
other observed galaxies and therefore that the gas responsible for the observed emission in the 119/im line resides 
predominantly in relatively dense and warm enviroments in comparison with these other sources. 

The other two observed emission lines, unlike the 119 /zm line, need not be coUisionally dominated. In the case of 
the H]^/2 3/2 — 1/2 163/im line, the most likely excitation mechanism is absorption of photons emitted by dust in the 
53 and 35 /im lines followed by radiative cascade. The upper level of this transition is 270 K above the ground state 
(Fig. EJ, so that excitation through collisions is expected to be ineffective in this line. The excitation mechanism of 
the Hi/2 — H3/2 1/2 — 3/2 79 /im line could be a mixure of collisional and radiative pumping. The upper level of this 
transition is 182 K above the ground state, so that a warm and dense region could, at least partially, excite the line 
through collisions. Nevertheless, the line could be also excited through the same infrared pumping mechanism that 
results in the observed 163 /im line emission. 

In conclusion, the 119 /im line is coUisionally excited, whereas absorption of photons emitted by dust in the 53 and 
35 /im lines probably dominates the excitation of the 163 /im line. The 79 /im OH line may in principle be excited 
through both mechanisms. 

5.3. Constraints on the spatial origin of the 119 fxm OH line 

In NGC 1068 two regions with very different physical conditions can account for the observed OH emission as 
discussed above: the compact nuclear region, and the ring and bar where the starburst is taking place. A warm and 
dense region is required to account for the observed 119 /im line emission, given that the line is coUisionally excited, so 
the warm and dense neutral nuclear region around the AG N should be considered a good candidate, despite its small 
size (~ 5"; e.g. iPlanesas et iiLlll99lHSchinnerer et alJl200nD . for the following reasons: 

(i) It is warm: there are ~ 10^ Mf:> of hot H2 (~ 20 00 K) distributed over ~ 5" l|Bhetz et al.lll994D that is thought 
to be UV or X-ray heated llRotaciuc et alJll99lD. Both PDRs and XDRs ca,n produce a range of temp eratures 
as high as a few times 10^ K (|Kaufman et al.lll999HSternberg fc Dalgarno1ll995t iMalonev et "311119961) . From 
the CO (4-3) to (1-0) line intensity ratio, iTacconi et al.l l)1994|) d erive ~ 80 K for the b u lk of the m olecular gas, 
with a mass of 3 x 10^ Mg enclosed in within the central 4" l|Helfer fc Blitz 111995(1 . iLutz et a l. (2000) have 
reported the detection of pure H2 rotational lines within the ISO-SWS aperture, and estimated ~ 2.5 x 10^ Mq 
at ~ 200 K, but these lines may also arise, at least partially, from the inner regions of the 3 kpc starburst ring. 

(ii) The molecular clouds within the nuclear region a re dense, although t here is some dispersion in the values 
derived by several authors based on HCN emission: ITacconi et al.l l|1994|) derived an H2 density of ~ 10'^ cm~^. 
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whereas subsequent observations and analysis bv lHelfer fc Blitz"! 1)1995(1 yielded a density of ~ 4 x 10® cm An 
int ermediate density o f ^ 5 x 10^ cm^'' from HON and CS, and lower for other tracers, has been recently derived 
bv lUsero et"al] l|200^ . 



(Hi) The OH abundance is expec ted to reach high values in regions exposed to strong inciden t UV field s (PDRs, 
[Sternberg & Dalearnol llQQfiD . and i n particular i n X-r ay dominated regions (XDRs, Lepp & Dalgarno 1996). 
The remarkable chemistry found bv Us ero et alJ 1)200 4) in the circumnuclear disk of NGC 1068 is indicative of 
an overall XDR and suggests a high OH abundance in the nuclear region. 

Given that the OH 119/im line is coUisionally excited, the possibility that the line might arise from the nuclear region 
can be checked by computing the amount of warm gas required to account for the observed emission: 



M^(Mo) = 1.6 X 10^ X 



10 



'5 1 



X(OH) 



5 X 10^ cm ^" 



4.3 X 10- 



< Q« > 



(1) 



n(H2) 

where X(OH) is the OH abundance relative to H2, and < ciu > is the collisional rate for excitation from the ground 
II3/2 3/2 level to the 113/2 5/2 one. Equation ^ assumes that, although the line could be optically thick, it is effectively 
optically thin, and makes use of the observed flux of 1.2 x 10~^^ erg s~^ cm~^. The reference value for the collisional 
rate, < q„ >= 4.3 x 10~^^ cm^ s~^, corresponds to gas at 80 K l)Offer et al. 111994^ : it decreases by a factor of « 2.7 
for gas at 50 K and increases by a factor of 3 for gas at 200 K. 

The reference OH abundance we use in this estimate, 10~^, is the result of two separate studies: first, calculations 
of molecular abundances by Lcdd & Dalearno (1996) have shown that the OH abundance in XDRs is expected to 
be about two orders of magnitude higher than the abundance of HCN and HCO^. The authors in fact suggested 
the possibility that the high HCN/CO ratio observed in the nuclear region of NGC 1068 could be a consequence 
of enhanced X- ray ionization. Sec ond, the possibility of a chemistry dominated by X-rays has found support from 
observations bv lTJsero et al.l 1)2004 ^ . who derive abundance ratios of HCN, HCO+ , and CN in general agreement with 
predictions for XDRs. Since the HCN abundance derived bv 'Usero et alJ l)2004D is ~ lO"'^, X(OH) in XDRs could 
attain values as high as 10~^. On the other hand, the density of 5 x 10^ cm~'^ derived bv lUsero et alJ l)2004|) has 
been adopted as the reference value in eq. ^ Finally, the mass of molecular ga s derived from the emission of several 
molecular tracers is expected to be at least ~ 2 x 10'' Mq ()Helfer fc Blitz lir995)) . which is similar to the value required 
in eq. From these estimates we conclude that, if the OH abundance is as high as ^ 10^^ (i.e. if the predictions 
for XDRs arc applicable to the nuclear region of NGC 1068), the bulk of the OH 119 /im line could arise there. This 
possibility would naturally explain why NGC 1068 is unique in its 119 fim line emission among galaxies with full LWS 
spectra. 

Finally we ask whether the OH 119/zm line could arise from an even more compact region, i.e., from a torus wit h 
a spatial scale of 1 pc surrounding the central AGN. According to typical parameters given bv lKrolik fc hevrii l)1989)l . 
a torus is expected to be hot (~ 10'^ K), could have densities of 10^ cm"'^, and therefore a mass of ~ 10^ M©. Also, 
the OH abundance is expected to be very high, 5 x 10~^ — 10"''. Eq. ^ shows that the relatively low mass of the 
torus (about 2 orders of magnitude lower than the entire nucleus) could be compensated by the the higher density, 
OH abundance, and temperature expected there, so that this possibility cannot be neglected. 

The reference values for the nuclear abundance and density adopted in eq. ^ £^re rather uncertain (and possibly 
extreme). The continuum models of §4 indicate that the mass associated with the 34 K dust component, which is 
identified with the starburst ring, is 1.2 x 10^ M0. If ~ 5% of this mass corresponds to warm molecular gas rich in OH, 
the amount of extended warm gas is 6 x 10^ Mq. According to eq. ^ the OH emission at 1 19 um can then also be 
expla i ned as arising in the ring if the associated PDRs, with assumed OH abundance of 2 x 10~® ()Sternberg fc Dalearno I 
[1^95"; 'Goicoechea fc Cernicharo ' '200^ iGonzalez- Alfonso et al. 1 1200^ . have densities of a few x 10^ cm~^. Since 
[PapadoDOulos fc Seaquist (1999b) found that most of the extended molecular gas resides in dense, compact clouds, 
this scenario seems also possible. However, the continuum from the starburst at 119 fim is strong, so that one expects 
that eq.^is in this case underestimating M^, and the quoted physical parameters, A'(OH) and ri(H2), are lower limits. 
The effect of dust emission is discussed in detail below. 

In conclusion, a definitive answer to the issue of the spatial origin of the OH 119 /im emission cannot be inferred 
from only the flux observed in the 119 /im line. Nevertheless, useful constraints on this subject are given: the line 
could be either explained as arising from the nucleus, with a required OH abundance ~ 10~^, or from the extended 
ring, with OH abundance > 2 x 10^® and density > a few x 10*^ cm^'^. Nevertheless, the radiative transfer models 
described below, which take into account the effect of the continuum emission and the excitation of the 79 and 163 
fim lines, point towards a nuclear origin of the OH emission. 

5.4. Outline of the models 

Analysis of the OH 79 and 163 /im lines requires the use of detailed radiative transfer calculations since, as pointed 
out above, the emission in these lines is expected to be strongly influenced by absorption of far-infrared continuum 
photons. We therefore proceeded to model the OH lines with two different codes, and confirmed that the results were 
in good agreement with each other. One of them, described in Gonzalez- Alfonso fc Cernicharo (1997, 1999), has been 
recently used to model the far-infrared spectrum of Arp 220 ( Gonzalez- Alfonso et al. II2004D . The other is a Monte 
Carlo radiative transfer code used as part of a detailed study of all the OH lines observed by ISO in galaxies l)Smith I 
ISmitheTalllmil . The code was developed for the Submillimeter Wave Astronomy Satellite (SWAS) mission 
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and is a modified and extended version of tlie Bernes code l|Bernes II1979I) but which includes d ust as weU as gas in 
the radiative transfer, and also corrects some optical depth calculations from the original code l)Ashbv et al.ll20 001 
Both methods are non-local, non-LTE, assume spherical symmetry, and include a treatment of continuum photons 
from dust mixed in with the gas. Also, both codes take input as a series of concentric shells, each of which is assigned 
a size, gas and dust temperature, H2 density, velocity and turbulent velocity width, and molecular abundance relative 
to H2 . The statistical equilibrium populations of OH in each spherical shell are computed by including the excitation 
by dust emission, excitation through collisions, and effects of line trapping. We ran two models to simulate the nucleus 
of the galaxy and the starburst extended ring, described in §5.4.1 and 5.4.2. 

5.4.1. Models for the nuclear emission: constraints on the spatial origin of the 79 and 163 iim OH lines 

We present models for the nuclear OH emission that implicitly assume that the 119 fim emission line arises from 
the nuclear region: X(OH)= 10~^ is adopted, as well as densities 5 x 10^ cm""^ for the bulk of the emitting gas. By 
assuming a pure nuclear origin for the 119 /im line, we check whether the other two OH lines could, in such a case, 
arise from the same nuclear region or require a more extended spatial origin. 

The models use the dust parameters derived for the nuclear molecular region described in §4. The possible contri- 
bution to the OH excitation of far-infrared photons arising from the components 1 & 2 of model AGN A is ignored, 
i.e. only the dust coexistent with the molecular gas is taken into account, with pre dicted flux of 2 4 — 27 Jy at 35-53 
/im (Fig. 0J). The gas temperature is assumed uniform and equal to 70 K (e.g. Tacconi et aLllTOQ^ . 

The H2 densities derived from the dust model (i.e. a peak density of < n/(H2) > = 500 cm""^ at the inner radius) 
are not compatible with the densities inferred from different molecular tracers. This indicates that t he medium 
is extremely clumped, as has been also argued elsewhere ('e.g. iCameron et aLlll993t iTacconi et al.lll994D . In order 
to account approximately for this dumpiness in our models, the following strategy is adopted: we use the "real" 
n(H2) ~ 5 X 10^ cm~^ values for the bulk of the gas, and compute the volume filling factor fy =< n(H2) > /n(H2), 
where the average value is that inferred from the dust model. The expected abundances of OH and the dust relative 
to H2, X(OH) and X(dust), are then multiplied by /„, so that the right OH and dust column densities are used in 
the calculations together with the right density values. The same density profile r~^ that was used in the dust model 
is adopted, so that /„ is uniform throughout the nuclear region. 

The modeled fluxes are convolved with the ISO-LWS grating resolution and are compared with the data in Fig. El 
Solid black lines show the results for the nuclear model that assumes Ar(OH)=10^^ and fy — 2x 10~^, the latter value 
implying a density in the outer regions (where the bulk of the emission is generated) of 5 x 10^ cm~^. The radial OH 
column density is iV(OH) = 2 x 10^^ cm~^. Besides the 119 /xm line, the model nearly reproduces the emission in the 
79 and 163 /im lines and is consistent with the upper limits given in Table El 

We also checked the excitation mechanism of the other OH lines by generating an additional model with the same 
parameters as above except for the continuum emission, which is now turned off. In this model, therefore, the lines 
are excited exclusively through collisions with H2. The resulting flux of the 119 /im line remains unchanged when 
the dust emission is ignored, confirming that the line is coUisionally excited. On the other hand, the flux densities of 
the 79 and 163 /xm lines decrease in the "pure-collisional" model by factors of 2 and 6, respectively, showing that the 
emission in these lines is much more affected by radiative pumping. We conclude that, if the OH abundance in the 
nucleus were high enough to account for the coUisionally excited 119 /im line, the observed fluxes in the 79 and 163 
/im lines can also be explained as arising in the same nuclear region. 

5.4.2. Models for the starburst emission 

Two simple different approaches have been used to model the OH emission from the starburst. First, we have 
roughly modeled the whole starburst region as a spherical shell with external radius of 1.5 kpc, thickness of 0.3 kpc, 
and average H2 density < n/(H2) > = 3. 7 cm ~'^, so that the corresponding continuum emission is reproduced with 
= 34 K (section0J. As shown in section l53l the OH 119 /im emission requires densities of a few x 10^ cm~^, so that 
we have assumed a volume filling factor f^ = 7.5 x 10^^ and therefore a "real" density n(H2) = 5 x 10^ cm~^. The 
kinetic temperature is assumed to be Tfc = 100 K, and the OH abundance is X(OH)= 2 x 10^^ x fh, where fh — 0.05 
is the assumed fraction of warm gas. The result of this model is shown in Fig. El (upper dotted lines). The 119 /im 
line is reproduced, but the flux densities of both the 79 /im and 163 /zm lines are strongly underestimated. The reason 
is that the model implicitly assumes that the continuum emission, responsible for the excitation of those lines, arises 
from a very large volume, so that the radiation density is weak and has negligible effect on the line excitation. 

Since the OH emission is expected to arise from compact, discrete PDRs in the vicinity of O or early B stars, 
where the continuum infrared radiation density is expected to be stronger than assumed above, we have also tried 
an alternative approach, which consists of modelling an individual "typical" cloud of the starburst. We first model 
the continuum from an individual cloud by assuming a central heating source and computing the dust equilibrium 
temperatures at each radial position that result from the balance of heating and cooling. The continuum model is 
adopted if (i) leaving aside a scaling factor (A^c, the number of clouds in the ensemble), the resulting SED is similar 
to that of the far-infrared emission of the starburst (i.e. the 34 K component found in sectional); the value of N,, is 
determined by requiring that the absolute continuum flux from the ensemble of clouds is equal to that observed for 
the 34 K component; and {ii) we require that the total mass of the ensemble does not exceed the mass inferred from 
the non-nuclear region (< 4 x 10^ Mq). Once the continuum is fitted, calculations for OH are performed by assuming 
Tfc = Td, and X(OH)= 2 x 10"^. 
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Several models with various density profiles were found to match the above two requirements. The common char- 
acteristic of all of them is the relatively high column density of the individual clouds, iV(H2) > lO^'^ cm~^, which is a 
consequence of the low effective dust temperature (34 K) of the SED. The results of the simplest model, characterized 
by a fiat density profile, n(H2) = 5 x 10^ cm~^, are given here for reference. With a radius of 5 x 10^^ cm, a stellar 
luminosity of 2 x 10^ L©, and A^^, = 4 x 10^, the resulting SED is similar to that of the 34 K component (§4). For 
these clouds we obtain a total mass of 1.8 x 10^ M©. The predicted OH emission/absorption is shown in Fig.El(lower 
dashed lines). In spite of the relatively high density and temperatures (28-250 K) throughout the cloud, the 119 fim 
line is predicted to be too weak, and the 79 fj,m line is predicted in absorption. We have found this result quite general: 
in models where the OH abundance is high enough and the radiation density becomes strong enough to pump the 163 
fim emission, the continuum at 79 and 119 /im is absorbed by OH and the predicted emission in the corresponding 
lines is reduced. Models that assume a density profile of generally predict the 119 /im line in absorption. In some 
models where the OH abundance was allowed to vary with radial position, the 79 /xm line was predicted in emission 
but by far too weak to account for the observed flux density. 

In conclusion, no starburst model is found to reproduce satisfactorily the emission observed in the three OH lines. If 
the local infrared radiation density is strong enough to pump the 163 /im line, the other two OH lines are expected to 
be weak or in absorption. Furthermore, the high density assumed for the starburs t region would p roduce a relatively 
high HCN/CO intensity ratio, which is on the contrary ^ 0.01 in the spiral arms ijHelfer fc Blitz l il995). Finally, the 
PDR models described in §3.2 indicate a density of 1 — 1.5 x 10^ cm~^, i.e. a density much lower than that required to 
account for the flux density of the OH 119 /im line. Therefore, and despite the simplicity of our models, taken together 
the analysis of the OH lines and the derived PDR parameters indicate that the bulk of the OH emission arises from 
X-ray dominated nuclear regions. 

6. CONCLUSIONS 
The main results of this article can be summarized as follows: 

• The complete far-infrared (50-200/im) spectrum of NGC 1068 has been observed for the first time. The far- 
infrared ISO-LWS spectrum has been complemented with the mid-infrared data of ISO-SWS and with shorter 
wavelength (UV, optical and near-IR) data from the literature to assemble a composite atomic spectrum as 
complete as possible with the aim of modeling the different line emission components at work. This approach 
has been necessary especially because of the poor spatial resolution of the ISO spectrometers, which were not able 
to spatially separate the emission components. The lines have been interpreted as arising from two physically 
distinct components: the AGN component and a starburst component, the first one nuclear and the second one 
located in the ring at a radius of 15-16 " from the nucleus. Both components are characterized by the presence 
of dust grains, producing strong continuum emission in the mid- and far-infrared. The density and ionization 
parameter of the ~ 5 x 10^ year old starburst are found to be uh ^ 100 cm~'^ and log U = -3.5, respectively. 
Three composite models h ave been computed with different AGN components: the first one has the ionizing 
continuum as derived from [Alexander et al.! (2000), showing a deep trough at energi es of a few Rydberg; the 
second has the monotonically decreasing ionizing continuum given bv lPier et al.l l)1994(l and the third has a "big 
blue bump". Two values of the visual extinction [Eb~v = 0.2 and 0.4) have been adopted to correct the optical 
and ultraviolet line fluxes for the reddening. The agreement betwe en the co mposite model with an AGN ionizing 
continuum characterized by the deep trough suggested by Alexand er et al.l |2000) is very satisfactory, taking into 
account both the simplicity of the photoionization models chosen to avoid dealing with too many free parameters 
and the large number of lines which originate in different physical regimes. The agreement between the observed 
spectrum and what is predicted using the canonical ionizing continuum is slightly poorer, while the presence of 
a big blue bump is ruled out. 

• The 50-200/^m continuum has been modeled using different components arising from both the nucleus and 
the starburst ring. For the nucleus, we have combined the dust emission from the ionized components in the 
narrow line regions modeled by CLOUDY with the neutral component reproduced by the radiative transfer code 
used for the OH molecular emission. For the starburst ring, our CLOUDY modelling of the ionized -I- PDR 
components could not reproduce the far- infrared emission, while similtaneously fitting the far-IR lines. Instead 
we fit the observed continuum by a neutral molecular component, reproduced by two gray body components at 
temperatures of 20K and 34K, assuming a steep (/3=2) dust emissivity law. 

• The unique OH emission in the 119 /im line cannot be explained in terms of OH excitation through absorption 
of 35 and 53 /xm photons emitted by dust, but rather it is coUisionally excited. This indicates the presence of a 
warm and dense region with high OH abundance. A simple excitation analysis yields two main alternatives for 
the spatial origin of the observed 119 /xm line emission: [i) the nuclear region, with 2 x 10 Mq of warm gas (80 
K), an average density of n(H2) = 5 x 10^ cm~'^, and an OH abundance of ~ 10~^; {ii) the starburst region, 
if ~ 5% of the associated mass (~ 6 x 10^ M©) is warm (~ 100 K), dense (a few x 10*^ cm~^), and rich in OH 
(A:(0H)- 2 X 10-^). 

• Radiative transfer models that simulate the emission/absorption in all the OH lines have been performed for 
both the nuclear and the starburst region. The models for the nucleus quantitatively account for the emission 
in the three OH lines if the nuclear physical conditions pointed out above are assumed. On the other hand, no 



10 



Spinoglio et al. 



starburst model is found to match the three OH hnes simultaneously, because the strong far-infrared continuum 
tends to produce absorption, or to weaken the emission, in the OH 119 and 79 /im lines (as observed in other 
galaxies). Therefore, although some contribution from the extended starburst cannot be ruled out, our models 
indicate that the bulk of the OH emission arises in the nuclear region. The high nuclear OH abundance required 
to explain the emission strongly suggest a chemistry deeply influenced by X-rays, i.e., an X-ray dominated region. 
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Fig. 1. — The AGN ionizing continua used as input for the photoionization model s of NGC 1068. The three continua differ in the 
frequency region betwe en l>Ej;,, j>100, while outside this region the Pier ct aQ 119941) spectrum was adopted. The solid line shows the 
continuum derived from lAlexander et al.l 120001) : the dashed line shows a simple power law interpolation; the dotted line shows the presence 
of the predicted "big blue bump". 
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Fig. 2. — The starburst spectral energy distributions used as input for the photoionization models of NGC 1068. The two continua are 
taken from Lcithcrcr e t al. ( 1993.) and represent a continuous starburst model (solid line) and an instantaneous model (broken line) , both 
with ages of 5 Myr. 
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Fig. 3. — The comparison of the composite models with the observations is shown as the ratio of modeled to observed flux ratio for each 
line, with the ionization potential in the x-axis. The assumed reddening is E(B-V)=0.2. Panels from top to bottom: model CMl, CM2 
and CMS. The short dashed lines represent flux ratios within a factor 3 either ways. 
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Fig. 4.— a) Spectral energy distribution of NGC 1068 and model fit. ISO-SWS fluxes are taken from Lutz et al. (200ff). The model fit 
is the composition of (i) the NLR components 1 & 2 of model AGN A, (ii) the emission from the molecular nuclear region, {Hi) the 34 K 
starburst emission, and (iv) the cold 20 K component, b) and c) Dust temperature versus the radial position and radial continuum opacity 
versus wavelength for the nuclear molecular region. 
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Fig. 5. — Energy level diagram of OH. Rotational levels with energies up to 600 K are shown; the three lines detected in NGC 1068 are 
indicated with arrows, as well as the 35 and 53 /xm lines that could play an important role in the radiative excitation. The wavelengths 
are indicated in //m. A-doubling is ignored because the A-doublets are not resolved with the ISO grating resolution. 



16 



Spinoglio et al. 



B 


7.5 


1 ^^^^ ' "3/2 5/2-3/2 ^ 














7 


^ Nucleus '''^^K^ 




o 




Starburst— ring 




6.5 


Starburst — clouds _ 

' 1 1 1 1 , 1 1 1 1 1 1 1 1 





116 118 120 122 




76 78 80 




160 162 164 166 




52 54 56 

X (yum) 



Fig. 6. — Comparison between the observed OH lines and model results. As indicated in the upper panel, the upper modeled spectrum 
(solid lines) corresponds to the model for the nucleus, the middle one (dotted lines) corresponds to the starburst modelled as a whole, and 
the lower one (dashed lines) corresponds to the starburst modelled as an ensemble of individual clouds (see text for details). 
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TABLE 1 

Measured line fluxes from the LWS and SWS grating spectra, with 1<t 

uncertainties. 



Line 


A 


Flux 


Aperture 


reference 




(urn) (10~ 


erg cm~'^) 


("') 






[Si IX] 3P2 ^ ^Pl 


2.584 


3.0 


14 


X 


20 


1 


[Mg VIII] 2p3/2 -> 2Pi/2 


3.028 


11. ± 1.1 


14 


X 


20 


1 


[Si IX] ^Pi^^Po 


3.936 


5.0 ± 0.6 


14 


X 


20 


1 


[Mg IV] 2Pi/2 ^ 2p3/2 


4.487 


7.6 ± 1.5 


14 


X 


20 


1 


[Ar VI] 2p3/2 ^ 2p^^2 


4.529 


15. ± 3. 


14 


X 


20 


1 


[Fe II] -"Pg/a ^ "''^9/2 


5.340 


5.0 


14 


X 


20 


1 


[Mg VII] 3P2 ^ 3Pl 
[Mg V] 3Pi ^ 


5.503 


13. 


14 


X 


on 


1 


5.610 


18. ± 2. 


14 


X 




1 


[Ar II] 2Pi/2 ^ 2pg^2 


6.985 


13. 


14 


X 


on 


1 


[Na III] 2p^^2 ^ ^-P3/2 


7.318 


5.8 


14 


X 


on 


1 


[Ne VI] 2p3/2^2p^^^ 


7.652 


110. ± 11. 


14 


X 


20 


1 


[Fe VII] ^ 


7.815 


3.0 


14 


X 


20 


1 


[Ar V] 3P2 ^ 3pj 


7.902 


< 12. 


14 


X 


20 


1 


[Na VI] 3p2 ^ 3p^ 


8.611 


< 16. 


14 


X 


20 


1 


[Ar III] 3Pi ^ 3P2 


8.991 


23.0 ± 3.3 


14 


X 


on 


1 


[Fe VII] 3P3 ^F2 


9.527 


4.0 


14 


X 


20 


1 


[S IV] ^P3/2 ^ 2p^^2 


10.510 


58. ± 6. 


14 


X 


20 


1 


[Ne II] 2P3/2 ^ 2Pi/2 


12.813 


70. 


14 


X 


07 


1 


[Ar V] 3pi ^ 3p|^ 


13.102 


< 16. 


14 


X 


07 


1 


[Ne V] 3P2 ^ 3pj 


14.322 


97. ± 9.7 


14 


X 


07 


1 


[Ne III] 3Pi ^ 3P2 


15.555 


160. ± 32. 


14 


X 


07 


1 


[Fe II] "*P7/2 - "^P9/2 


17.936 


< 10. 


14 


X 


07 


1 


[S III] 3P2 ^ 3p^ 


18.713 


40. 


14 


X 


07 
2 / 


1 


[Ne V] 3pi ^ 3p|j 


24.317 


70. ± 7. 


14 


X 


27 


1 


[O IV] 2P3/2 ^ 2p^^2 


25.890 


190. ± 20. 


14 


X 


27 


1 


[Fe II] "6D7/2 ^ "''Og/z 


25.988 


8. 


14 


X 


27 


1 


[S III] ^'Pi ^ ^Po 


33.481 


55. 


20 


X 


33 


1 


[Si II] 2P3/2 ^ ^P^,2 


34.814 


91. 


20 


X 


33 


1 


[Ne III] 3Po ^ ^Pi 


36.013 


18. 


20 


X 


33 


1 


[O III] 3p2 ^ 3pj 


51.814 


114. ± 3. 




80 




2 


[N III] 2P3/2 ^ ^P^/2 


57.317 


51.4 ± 2.5 




80 




2 


[O I] 3Pi ^ 3p2 


63.184 


156. ± 1. 




80 




2 


[O III] 3Pi ^ 3Po 


88.356 


111. ± 1. 




80 




2 


[N II] 3P2 ^ 3p^ 


121.897 


30.5 ± 1.1 




80 




2 


[O I] 3Po ^ 3p^ 


145.525 


11.9 ± 0.4 




80 




2 


[C II] 2p3/2 ^ 2p^^2 


157.741 


216. ± 1. 




80 




2 


OH 2ni/25/2-2n3/23/2 

OH 2ni/2i/2-2ni/23/2 

OH 2n3/25/2-2n3/23/2 


34.60/34.63 


< 3. 


20 


X 


33 


2 


79.11/79.18 


14.4 ± 1.5 




80 




2 


119.23/119.44 


11.9 ± 1.2 




80 




2 


OH 2ni/23/2-2ni/2i/2 


163.12/163.40 


7.42 ± 0.65 




80 




2 


Note. — fll: fromlLutz et 


alJ I2000H and, where 


errors arc available. lAlexander et 


alJ 120001): 



(2) : this work 
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TABLE 2 

Comparison of observed line fluxes with AGN model predictions 





Line id. A 






Flux (10^13 


erg s ^ cm 










(Min ) 


Observed/D^ /D^ 


AGN A model» 


AGN B model4 


AGN C models 








Comp. 1 


Comp. 2 


Comp. 1 


Comp. 2 


Comp. 1 


Comp. 2 


O VI A .1032+. 1037 


*H7 4/4'HS4 /402 


32.6+19.8 


2.28+1.94 


8.26+5.36 


0.31+0.27 


56.7+31.2 


3.41+2.10 


(Lyo) 


n A .1215 


101. 8/3562. /602 


179 


619 


149. 


500 


239. 


638. 


N IVl A .1487 


5 1 /1 03 /22 9 


25.2 


6.64 


36.3 


6.28 


59.6 


47.5 


(CIV) 


n A .1549 


39 7/790 /1 77 


142. 


77.6 


110.9 


56.5 


286. 


329. 


Hell A .1640 


21 4/426 /q5 5 


112. 


85.9 


108 


81.8 


346. 


345. 


[Ne V] A .3426 


15. 7/95. /38. 7 


97.7 


23.3 


68.2 


7.36 


271. 


115 


[Ne III] A .3869+. 3968 


19 2/97 /43 2 


37.1+11.2 


41.5+12.5 


102.8+31.0 


102. +30.7 


81.6+24.6 


143. +43.0 


Hell A .4686 


fi 1 /27 6/13 


15.0 


12.3 


14.7 


11.8 


46.5 


48.4 




O III] A .4959+. 5007 


256 /Qfi4 /496 

^'jyj. / cut:. / '-t.jyj 


86.9+262 


99.5+294 


197.+593 


185. +557. 


148.+446. 


241+726 




Si VI] A 1.96 


8 O/Q 2/8 6 


11.2 


1.53 


10.8 


1.86 


23.3 


3.74 




Si VII] A 2.48 


8.3 


6.41 


0.1 


9.12 


0.18 


11.0 


0.22 




Si IX 


A 2.584 


3.0 


0.49 


— 


0.84 


— 


0.32 


— 




Mg VIII] A 3.028 


11. 


2.97 


— 


4.80 


— 


4.29 


— 




Si IX 


A 3.936 


5.4 


0.89 


— 


1.56 


— 


0.56 


— 




Mg IV] A 4.487 


7.6 


3.85 


7.98 


4.29 


6.92 


9.80 


19.5 




Ar VI] A 4.529 


15. 


10.3 


4.0 


13.4 


3.10 


26.2 


20.1 




Mg VII] A 5.503 


13. 


9.68 


0.13 


11.4 


.075 


20.0 


0.45 




Mg V] A 5.610 


18. 


10.7 


8.63 


7.85 


3.84 


34.6 


29.5 




Ar 11] 


A 6.985 


13. 


4.56 


7.53 


0.93 


1.36 


0.87 


1.63 




Na III] A 7.318 


5.8 


0.56 


0.59 


1.27 


1.51 


0.67 


1.21 




Ne VI] A 7.652 


110. 


153.7 


11.9 


105.4 


3.2 


393. 


62.4 




Fe VII] A 7.815 


3.0 


1.83 


1.63 


1.65 


0.71 


6.34 


6.30 




Ar V] A 7.902 


< 12. 


2.14 


1.83 


2.62 


1.92 


3.72 


6.52 




Na VI] A 8.611 


< 16. 


1.20 


0.16 


0.95 


— 


2.95 


0.65 




Ar III] + [Mg VII] A 8.991 


25. 


4.63+11.7 


7.34+0.17 


5.97+14.1 


12.4+0.10 


4.83+23.8 


10.7+0.58 




Fe VII] A 9.527 


4.0 


7.50 


7.0 


6.75 


3.08 


25.8 


26.6 




S IV] A 10.510 


58. 


38.4 


60.5 


85.3 


126. 


75.8 


187. 




Ne II 


A 12.813 


70. 


5.86 


20.8 


2.20 


4.44 


1.62 


3.26 




Ar V 


A 13.102 


< 16. 


2.23 


2.73 


2.75 


2.88 


3.93 


9.57 




Ne V 


A 14.322 


97. 


91.4 


66.8 


83.2 


28.7 


270. 


273. 




Ne III] A 15.555 


160. 


44.4 


52.0 


76. 


110. 


42.0 


86.9 




S III] A 18.713 


40. 


25.8 


56.0 


21.7 


75.7 


18.4 


64.5 




Ne V] A 24.317 


70. 


40.8 


52.6 


35.8 


22.1 


116. 


215. 




O IV] A 25.890 


190. 


24.1 


71.9 


28.3 


69.7 


56.4 


240. 




S III] A 33.481 


55. 


8.00 


29.8 


5.57 


41.1 


4.84 


38.4 




Si II] A 34.814 


91. 


12.9 


22.6 


6.01 


24.2 


5.28 


29.4 




Ne III] A 36.013 


18. 


3.32 


4.4 


5.74 


9.40 


3.20 


7.46 




O HI 


A 51.814 


110. 


9.04 


33.9 


16.3 


63.2 


9.06 


48.5 




N ni] A 57.317 


51. 


2.23 


11.4 


2.83 


18.9 


1.88 


16.1 




O I] A 63.184 


156. 


3.84 


1.8 


1.57 


2.26 


1.24 


3.10 




O HI 


A 88.356 


110. 


1.22 


8.2 


2.25 


15.5 


1.23 


12.6 




N II] A 121.897 


30. 


0.28 


0.87 


.07 


0.40 


.06 


0.52 




b I] A 145.525 


12. 


0.24 


0.14 


.09 


0.18 


.07 


0.24 




C II] A 157.741 


220. 


0.50 


1.58 


0.18 


0.97 


.14 


1.09 



Dercddcncd line flux, assuming E^ — v — 0.4 
■^DcrcddcnGd line flux, assuming Eb-v — 0.2 

^AGN A param eters: component 1: Log U— -1., Log n— 4, internal radius c::^ 21 pc, external radius Ci^ 109 pc, ionizing spectrum from 
lAlexander et al ] 120011): com ponent 2: Log U=-2., Log n— 3.3, internal radius 2:; 153 pc, external radius 2± 362 pc, ionizing spectrum 
from lAlexander^Trr i2000y 
'^AGN B parameters: same as AGN A models, but with the ionizing spectrum from lPier et all il994i) 
^AGN C parameters: same as AGN A models, but with the ionizing spectrum that includes a big blue bump (see text) 
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TABLE 3 

Comparison of observed line fluxes with the Ring Starburst model predictions 



Line id. A Flux (10 erg s ^ cm ^) 



(/im ) 


Observed 


SBR Ai 


SBR B2 


SBR & 


SBR 


SBR EB 


SBR pS 


O VI A .1032+. 1037 


37.4/4334./402. 


— 


— 


— 


— 


— 


— 


(Lya)n A .1215 


101.8/3562./602 


1330. 


865. 


997. 


875. 


1010. 


917. 


N IV] A .1487 


5.1/103./22.9 


71.9 


0.20 


85.8 


.020 


93.7 


0.03 


(ClV)n A .1549 


39.7/790./177. 


795. 


0.70 


858. 


0.80 


924. 


0.87 


HcII A .1640 


21.4/426./95.5 


693. 


50.2 


696. 


50.5 


696. 


50.8 


Nc V] A .3426 


15.7/95./38.7 


58.1 


— 


69.6 


— 


73.9 


— 


Nc III] A .3869+. 3968 


19.2/97./43.2 


314.+94.7 


31.1+9.4 


309. +93.1 


31.5+9.5 


316. +95.7 


32.+9.7 


HcII A .4686 


6.1/27.6/13. 


97.3 


7.4 


97.7 


7.4 


97.3 


7.4 


O III] A .4959+.5007 


256./964./496 


1330.+3990. 


25.1+75.6 


1383+4158. 


26.0+78.2 


1412. +4257. 


26.8+80. 


Si VI] A 1.96 


8.0/9.2/8.6 


— 


— 


— 


— 


— 


— 


Si VII] A 2.48 


8.3 


— 


— 


— 


— 


— 


— 


Si IX] A 2.584 


3.0 


— 


— 


— 


— 


— 


— 


Mg VIII] A 3.028 


11. 


— 


— 


— 


— 


— 


— 


Si IX] A 3.936 


5.4 


— 


— 


— 


— 


— 


— 


Mg IV] A 4.487 


7.6 


3.30 


— 


3.50 


— 


3.53 


— 


Ar VI] A 4.529 


15. 


6.20 


— 


7.72 


— 


8.02 


— 


Mg VII] A 5.503 


13. 


— 


— 


— 


— 


— 


— 


Mg V] A 5.610 


18. 


3.10 


— 


3.63 


— 


3.73 


— 


Ar II] A 6.985 


13. 


5.90 


3.70 


4.29 


3.70 


4.12 


3.73 


Na III] A 7.318 


5.8 


0.60 


0.10 


0.54 


0.11 


0.54 


0.11 


Ne VI] A 7.652 


110. 


16.9 


— 


21.1 


— 


21.8 


— 


Fe VII] A 7.815 


3.0 


0.40 


— 


0.48 


— 


0.50 


— 


Ar V] A 7.902 


<12. 


7.50 


— 


8.84 


— 


9.11 


— 


Na VI] A 8.611 


<16. 


— 


— 


— 


— 


— 


— 


Ar III] + [Mg VII] A 8.991 


25. 


46.5 


7.10 


43.2 


7.13 


42.2 


7.13 


Fc VII] A 9.527 


4.0 


1.70 


— 


2.04 


— 


2.11 


— 


S IV] A 10.510 


58. 


312. 


0.80 


346. 


0.84 


333. 


0.81 


Nc II] A 12.813 


70. 


7.80 


5.50 


5.44 


5.38 


5.38 


5.61 


Ar V] A 13.102 


<16. 


21.9 


— 


25.5 


— 


24.9 


— 


Nc V] A 14.322 


97. 


271. 


— 


315. 


— 


314. 


— 


;Nc III] A 15.555 


160. 


295. 


50.5 


276. 


50.2 


270. 


49.8 


S III] A 18.713 


40. 


180. 


19.2 


178. 


20.3 


203. 


23.7 


Nc V] A 24.317 


70. 


303. 




346. 




294. 




OIV] A 25.890 


190. 


1660. 


1.90 


1762. 


2.00 


1419. 


1.62 


S 111] A 33.481 


55. 


785. 


87.5 


673. 


81.2 


350. 


46.2 


Si 11] A 34.814 


91. 


126. 


54.8 


89.4 


49.5 


41.9 


29.4 


Ne III] A 36.013 


18. 


52.8 


8.90 


49.2 


8.91 


47.5 


8.74 


O III] A 51.814 


110. 


2980. 


81.5 


3217 


87.5 


2524. 


71.9 


N III] A 57.317 


51. 


1110. 


40.0 


921 


35.0 


372.9 


14.7 


O I] A 63.184 


156. 


165. 


130. 


113. 


127. 


116. 


124. 


O III] A 88.356 


110. 


4980. 


138. 


3729. 


109. 


977. 


29.0 


N II] A 121.897 


30. 


80.0 


37.3 


47.8 


29.7 


12.3 


8.78 


O I] A 145.525 


12. 


16.2 


12.9 


11.0 


12.5 


10.1 


11.3 


C II] A 157.741 


220. 


739. 


465. 


181. 


172. 


67.6 


46.5 



Note. — *: this line wsis used for normalization 
"'"SBR A parameters: Log U— -2.5, Log n— 1.0, ionizing speetrum from Starburst99 with instantaneous star-formation law, M — IO^Mq, 
IMF: —2.35 M^p — IOOMq, Miow — '^^q, nebular emission included, Z— 0.020, age of 5 Myr. The integration was stopped at a 
temperature of 50K, the adopted abundances are those relative to HII regions and grain emission is included. The adopted number of 
clouds is 33000. 

■^SBR B parameters: Log U=-3.5, Log n=1.0, all other parameters as for SBR A. The number of clouds adopted is 33000 
^SBR C parameters: Log U— -2.5, Log n— 2.0, all other parameters as for SBR A. The number of clouds adopted is 3300 
^SBR D parameters: Log U— -3.5, Log n— 2.0, all other parameters as for SBR A. The number of clouds adopted is 3300 
^SBR E parameters: Log U— -2.5, Log n— 3.0, all other parameters as for SBR A. The number of clouds adopted is 330 
^SBR F parameters: Log U=-3.5, Log n=3.0, all other parameters as for SBR A. The number of clouds adopted is 330 
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TABLE 4 

Comparison of observed line fluxes with composite model 

rnEnicTioxH 



Line id. A 


r lux l^iU 


13 -1 
erg s 


cm J 




(/im ) 


Observed/D /D 


CMl^ 


CM2 


CM3 


O VI A .1032+.1037 


37.4/4334. /402. 


56.6 


14.2 


93.4 


(Lya)n A .1215 


101. 8/3562. /602 


1673. 


1524. 


1752. 


N IV] A .1487 


5. 1/103. /22. 9 


31.8 


42.6 


107. 


(ClV)n A .1549 


39. 7/790. /177. 


220. 


168. 


616. 


HcII A .1640 


21. 4/426. /95. 5 


248. 


240. 


741. 


[Nc V] A .3426 


15. 7/95. /38. 7 


121. 


75.6 


386. 


[Nc III] A .3869+.3968 


19.2/97./43.2 


143. 


307. 


333. 


HcII A .4686 


6 1 /27 6/13 


34.7 


33.9 


102.3 




O III] A .4959+.5007 


256 /qfi4 /4Q6 


847. 


1636. 


1665. 




Si VI] A 1.96 


8 0/9 2/8 6 


12.7 


12.7 


27.0 




Si VII] A 2.48 


8.3 


6.5 


9.3 


11.2 




Si IX] A 2.584 


3.0 


0.5 


0.8 


0.3 




Mg VIII] A 3.028 


11. 


3.0 


4.8 


4.3 




Si IX] A 3.936 


5.4 


0.9 


1.6 


0.6 




Mg IV] A 4.487 


7.6 


11.8 


11.2 


29.3 




Ar VI] A 4.529 


15. 


14.3 


16.5 


46.3 




Mg VII] A 5.503 


13. 


9.8 


11.5 


20.5 




Mg V] A 5.610 


18. 


19.3 


11.7 


64.1 




Ar II] A 6.985 


13. 


15.8 


6.0 


6.2 




Na III] A 7.318 


5.8 


1.3 


2.8 


2.0 




INc VIJ A 7.652 


110. 


166. 


109. 


455. 




Fe VII] A 7.815 


3.0 


3.5 


2.4 


12.6 




Ar V] A 7.902 


<12. 


4.0 


4.5 


10.2 




Na VI] A 8.611 


<16. 


1.4 


0.9 


3.6 




Ar III] + [Mg VII] A 8.991 


25. 


31. 


39.7 


47.0 




Fc VII] A 9.527 


4.0 


14.5 


9.8 


52.4 




S IV] A 10.510 


58. 


99.7 


212. 


263. 




Nc II] A 12.813 


70. 


32.0 


12.0 


10.3 




Ar V] A 13.102 


<16. 


5.0 


5.6 


13.5 




Nc V] A 14.322 


97. 


158. 


112. 


543. 




Nc III] A 15.555 


160. 


147. 


236. 


179. 




S III] A 18.713 


40. 


102. 


118. 


103. 




Nc V] A 24.317 


70. 


93.4 


57.9 


331. 




O IV] A 25.890 


190. 


98. 


100. 


298. 




S III] A 33.481 


55. 


119. 


128. 


124. 




Si II] A 34.814 


91. 


85. 


79.7 


84.2 




Ne III] A 36.013 


18. 


16.6 


24.0 


19.6 




O III] A 51.814 


110. 


130. 


167. 


145. 




N III] A 57.317 


51. 


48.6 


56.7 


53. 




O I] A 63.184 


156. 


133. 


131. 


131. 




O III] A 88.356 


110. 


118. 


127. 


123. 




N II] A 121.897 


30. 


30.8 


30. 


30.3 




O I] A 145.525 


12. 


12.9 


12.8 


12.8 




C II] A 157.741 


220. 


174. 


173. 


173. 



Dereddened line flux, assuming Eb-v = 0.4 
■^Dereddened line flux, assuming Es-v = 0.2 
^CMl = AGN A + SBR D 
= AGN B + SBR D 
5CM3 = AGN C + SBR D 



TABLE 5 

Comparison of observed line fluxes with model 
predictions for the nuclear region 



Line id. A 


Flux (10-12 
Observed 


erg s ^ cm 
Modeled 


Notes 


34/im 


< 0.1 


-0.43 


(absorption) 


53/im 


< 1.2 


-0.66 


(absorption) 


79//m 


1.1 


1.13 




84//m 


< 1.2 


0.06 




98/[fm 


< 1.2 


0.17 




119/um 


1.3 


1.60 




163/um 


0.38 


0.35 





